% syms a1 a2 a3 a4
% a=[a1 a2 0 0; a3 a4 0 0; 0 0 1 0; 0 0 0 1];
% b=sym(toym(4));
% b

a1=diag([1 2 3]);
a2=diag([1 1/2 1/3]);
a3=[0 0 1; 0 2 0 ; 3 0 0];
a4=[0 0 1; 0 1/2 0 ; 1/3 0 0];
b=toym(3);
